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Abstract 

We study analytically the asymptotic late-time evolution of realistic rotat- 
ing collapse. This is done by considering the asymptotic late-time solutions 
of Teukolsky's master equation, which governs the evolution of gravitational, 
electromagnetic, neutrino and scalar perturbations fields on Kerr spacetimes. 
In accordance with the no-hair conjecture for rotating black-holes we show 
that the asymptotic solutions develop inverse power-law tails at the asymp- 
totic regions of timelike infinity, null infinity and along the black-hole outer 
horizon (where the power-law behaviour is multiplied by an oscillatory term 
caused by the dragging of reference frames) . The damping exponents charac- 
terizing the asymptotic solutions at timelike infinity and along the black-hole 
outer horizon are independent of the spin parameter of the fields. However, the 
damping exponents at future null infinity are spin dependent. The late-time 
tails at all the three asymptotic regions are spatially dependent on the spin 
parameter of the field. The rotational dragging of reference frames, caused by 
the rotation of the black-hole (or star) leads to an active coupling of different 
multipoles. 
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I. INTRODUCTION 



The no-hair theorem, introduced by Wheeler in the early 1970s [|l], states that the ex- 
ternal field of a black-hole relaxes to a Kerr-Newman field characterized solely by the black- 
hole's mass, charge and angular-momentum. The mechanism for this relaxation process of 
neutral fields was first analyzed by Price M for a nearly spherical collapse. The physical 
mechanism by which a charged black-hole, which is formed during a gravitational collapse 
of a charged matter, dynamically sheds its charged hair was first studied in P,ffl. However, 
these analysis were restricted to spherically symmetric backgrounds, i.e., to the Schwarschild 
and Reissner- Nordstrom black-holes. On the other hand, the physical process of stellar core 
collapse to form a black-hole is expected to be non-spherical in nature because of stellar ro- 
tation. The analogous problem of the dynamics of massless perturbations outside realistic, 
rotating black-holes is much more complicated due to the lack of spherical symmetry. 

Recently, this problem was addressed numerically and analytically for the case of 
scalar-perturbations 0. However, the evolution of higher-spin fields (and in particular grav- 
itational perturbations) outside realistic, rotating black-holes was not analyzed analytically 
so far. In addition, the works done so far were restricted to the asymptotic regions of 
timelike infinity and along the black-hole outer- horizon |7]]. The asymptotics at null 
infinity of rotating collapse was not studied so far. In this paper we give for the first time 
a full analytic analysis of the late-time evolution of realistic rotating collapse. Our analysis 
considers massless fields with arbitrary spin (and in particular gravitational perturbations). 
In addition, we consider the late-time evolution at all the three asymptotic regions: timelike 
infinity, future null infinity and along the black-hole outer horizon. 

There are two different approaches to the study of perturbations of Kerr spacetimes. 
The first is to consider perturbations of the metric functions. However, this direct approach 
leads to gauge-dependent formulations. An alternative approach is to consider curvature 
perturbations (perturbations of the Weyl scalars). Based on the tetrad formalism by New- 
man and Penrose, Teukolsky derived a master equation governing the perturbations of Kerr 
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spacetimes |8]|§. The main goal of this paper is to analyze the asymptotic late-time evolution 
of a realistic rotating collapse. This is done by studying the asymptotic late-time solutions 
of the Teukolsky equation. 

The plan of the paper is as follows. In Sec. || we give a short description of the phys- 
ical system and formulate the evolution equation considered. In Sec. fTJ we formulate the 
problem in terms of the black-hole Green's function using the technique of spectral decom- 



position. In Sec. [TV] we study the late-time evolution of perturbations of Kerr spacetimes, 
i.e., the asymptotic late-time solutions of the Teukolsky equation. We conclude in Sec. [V| 
with a brief summary of our results. 



II. DESCRIPTION OF THE SYSTEM 



We consider the evolution of gravitational, electromagnetic, neutrino and scalar (mass- 
less) perturbations fields outside a rotating collapsing star. The external gravitational field 
of a rotating object of mass M and angular-momentum per unit-mass a is given by the Kerr 
metric, which in Boyer-Lindquist coordinates takes the form 
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where E = r 2 + a 2 cos 2 9 and A = r 2 — 2Mr + a 2 . Throughout this paper we use G — c — 1. 

Using the Newman-Penrose formalism, Teukolsky §|| derived a master equation that 
governs the evolution of perturbations of the Kerr spacetime 
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where the parameter s is the spin-weight of the field. The Teukolsky's master equation 
is valid for scalar perturbations fields (s = 0), neutrino perturbations fields (s = ±1/2), 



electromagnetic perturbations fields (s = ±1) and gravitational perturbations (s = ±2). 
The field quantities ip which satisfy this equation (for the various values of the spin parameter 
s) are given in M. Resolving the field in the form 
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one obtains a wave-equation for each value of m 
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where the tortoise radial coordinate y is defined by dy = I —j^-dr. The coefficients Bi(r,6) 
are given by 
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III. FORMALISM 



The time-evolution of a wave-field described by Eq. (^) is given by 

t) = 2n [ f {B 1 {z') \G{z, z'; t)%{z', 0) + G t (z, z'; 0) 



B 2 {z')G{z, z'\ O^sinO'de'dy' 



(8) 



for t > 0, where z stands for (y, 6 1 ). The (retarded) Green's function z'\ i) is defined by 
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The causality condition gives us the initial condition G(z, z'\ t) = for t < 0. In order to 
find G(z, z'\ t) we use the Fourier transform 
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where s S' z m (6 l ,aw) are the spin- weighted spheroidal harmonics which are solutions to the 
angular-equation || 
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For the aw = case, the eigenfunctions s S™(0,aw) reduce to the spin-weighted spherical 
harmonics s Y{ n {9, 0) = s SJ n (9)e tm ' p , and the separation constants s A 7 l l {aw) are simply gA™ = 
(l-s)(l + s + l) 0. 

The Fourier transform is analytic in the upper half u>-plane and it satisfies the equation 
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where K = (r 2 +a 2 )w—am, X = A+a 2 w 2 —2amw and H = s(r— M)/(r 2 +a 2 )+rA/(r 2 +a 2 ) 2 . 
G(z, z'; t) itself is given by the inversion formula 
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where c is some positive constant and Iq = max(\m\, \s\). 

Next, we define two auxiliary functions *&i(z, w) and ^2{z,w) which are (linearly inde- 
pendent) solutions to the homogeneous equation 
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The two basic solutions that are required in order to build the black-hole Green's function 
are defined by their asymptotic behaviour: 
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where k = w — mw + , w + = a/(2Mr + ). Let the Wronskian be 
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where W(w) is y-independent. Thus, using the two solutions and \^2, the black-hole 
Green's function can be expressed as 
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In order to calculate G(z,z';t) using Eq. ([13]), one may close the contour of integra- 
tion into the lower half of the complex frequency plane. Then, one finds three distinct 
contributions to G(z, z'\ t) 



1. Prompt contribution. This comes from the integral along the large semi-circle. This 
term contributes to the short-time response. 

2. Quasinormal modes. This comes from the distinct singularities of G(y,y';w) in the 
lower half of the complex w-plane, and it is just the sum of the residues at the poles of 
G(y, y'\ w). Since each mode has Imw < this term decays exponentially with time. 

3. Tail contribution. The late-time tail is associated with the existence of a branch cut 



(in \l/ 2 ) ]TT] , usually placed along the negative imaginary w-axis. This tail arises from 
the integral of G(y,y';w) around the branch cut and is denoted by G c . As will be 
shown, the contribution G c leads to an inverse power-law behaviour (multiplied by a 
periodic term along the black-hole outer horizon) of the field. Thus, G c dominates 
the late-time behaviour of the field. 



The present paper focuses on the late-time asymptotic solutions of the Teukolsky equa- 
tion. Thus, the purpose of this paper is to evaluate G c (z, z'\ t). 



IV. LATE-TIME BEHAVIOUR OF REALISTIC ROTATING COLLAPSE 



A. The large-r (low-w) approximation 



It is well known that the late-time behaviour of massless perturbations fields is deter- 
mined by the backscattering from asymptotically far regions [12,!]]. Thus, the late-time 
behaviour is dominated by the /oiofrequencies contribution to the Green's function, for only 
low frequencies will be backscattered by the small effective potential (for r 3> M) in Eq. 
(14). Thus, as long as the observer is situated far away from the black-hole and the initial 
data has a considerable support only far away from the black-hole, a large-r (or equivalently, 
a low-w) approximation is sufficient in order to study the asymptotic late-time behaviour of 
the field jl3 |. Expanding Eq. (]TJ]) for large r one obtains [using A = (/ — s)(Z + s + 1) +0(aw) 
and neglecting terms of order O(^) and higher] 
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It should be noted that this equation is the correct form of equation (66) of Ref. 
term s 2 /r 2 should not appear in that equation). 

We introduce a second auxiliary field defined by 
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where x = —2iwr. <p{x) satisfies the confluent hypergeometric equation 
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Thus, the two basic solutions required in order to build the Green's function are (for r ^> 
MM) 



4f! = Ar l+1 e iwr M(l + s + 1 - 2iwM, 21 + 2, -2iwr) , (23) 



and 



# 2 = Br i+1 e lwr U(l + 8 + 1- 2iwM, 21 + 2, -2iwr) , (24) 

where A and B are normalization constants. M(a, b, z) and U(a, b, z) are the two standard 
solutions to the confluent hypergeometric equation []14| |. From these solutions it is clear that 



the black-hole's rotation parameter a is irrelevant in the context of the asymptotic form of 
the Green's function (to leading order in the inverse time). However, as will be shown below, 
the rotational dragging of reference frames, caused by the rotation of the black- hole (or star), 
has (two) important effects (which are absent in the non-rotating case) on the asymptotic 
late-time evolution of rotating collapse. Hence, in order to find the asymptotic form of the 
Green's function it is sufficient to analyze the late-time solutions in the spherically symmetric 
limit (i.e, on the Schwarchild background). 

A note is needed here on the relations between the various equations governing the 
perturbations of the Schwarchild spacetime. There are two different approaches to the 
study of linearized gravitational perturbations: by considering metric perturbations, or, 
alternatively, by considering perturbations of the Weyl scalars. Equations governing metric 
perturbations of the Schwarchild spacetime were derived by Regge and Wheeler |15| (for odd- 
parity perturbations) and by Zerilli (for polar perturbations). An alternative approach, 
based on the tetrad formalism by Newman and Penrose was used by Teukolsky [7,8] to derive 
a separable wave equation for perturbations of the Weyl Scalars. These scalars, constructed 
from the Weyl tensor in a given tetrad basis, characterize the gravitational field in vacuum. 
They allow a more convenient approach to the study of gravitational perturbations, due 
to their scalar nature. As was first shown by Chandrasekhar flT7| ], the a — > limit of the 



Teukolsky equation, the Bardeen-Press equation p8| , can easily be transformed into the 



Regge- Wheeler equation. However, this equivalence is restricted to the frequency domain. 
In other wards, any solution (in the frequency domain) of the Bardeen-Press equation can 
be transformed in a trivial manner into a solution of the Regge- Wheeler equation, and vise 
versa. However, these solutions have, of coarse, a different w- dependence (as can be verified 
for example from their asymptotic forms Eqs. ( P3"D and ([Zip). Hence, this difference in the 
w-dependence of the solutions, when integrated in the complex frequency plane (in order to 
obtain the temporal dependence of the solutions) may lead to a non-trivial difference in the 
time dependence of the solutions. Thus, it is not apriori guaranteed that the asymptotic 
late-time solutions of Teukolsky's equation and the Regge- Wheeler equation (studied by 
price for the spherically symmetric case) will have the same temporal behaviour (Indeed, 
it is found that the damping exponents, describing the fall-off of the asymptotic solutions 
at future null infinity are spin-dependent, a phenomena not found for the Regge- Wheeler 
asymptotic solutions, see Sec. [IV C| below). Furthermore, one can easily verify that even 
the a —>■ limit of the Teukolsky equation depends on the spin parameter s [see Eq. (0) 
and its asymptotic form Eq. QI9"D]. Hence, it is not surprising that the asymptotic late-time 
solutions of the Teukolsky equation will turn out to be spin-dependent. 

The function U(a,b,z) is a many-valued function, i.e., there will be a cut in \t> 2 - Using 
Eq. fll3D , one finds that the branch cut contribution to the Green's function is given by 

1 00 p—ioo _ 
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s Si{8,aw) s S l {8',aw)e- iwt dw . (25) 

(For simplicity we assume that the initial data has a considerable support only for r-values 
which are smaller than the observer's location. This, of course, does not change the asymp- 
totic Za£e-tiine behaviour). 

Using the fact that M(a,b,z) is a single- valued function and Eq. 13.1.10 of |L4]] (taking 
the b — > k limit, where k is an integer), one finds 

fi(r,TO 2 ")=fi(r,w), (26) 
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Since iy(u>) is r-independent, we may use the large-r asymptotic expansions of the confluent 
hypergeometric functions (given by Eqs. 13.5.1 and 13.5.2 in [Q) in order to evaluate it. 
One finds 

W ^=* T { 1 + s\i-Jm)2^ - (30) 
(Of coarse, using the \z\ — > limit of the confluent hypergeometric functions, we obtain the 



same result). Thus, substituting Eqs. (p9|) and (|30|) in Eq. (|25|) one finds 
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As was explained, the late-time behaviour of the field should follow from the Zoiy-frequency 
contribution to the Green's function. It is clear that all the rot at ion- dependent terms of Eq. 
(PI) can be neglected in the aw — > limit. Hence, one may replace the functions s S^ !l {9 1 aw) 
by the spin- weighted spherical harmonics s Y{ n {9,(j)). In addition, we use the approximation 

which is valid for w — > 0. With these substitutions Eq. ( J3T]) becomes 

{ ' , ] ~h vrA 2 [(2/ + l)!] 2 sYl{9)sYl {d) 

/ ^ 1 (y,w)^ 1 (y',w)w 2l+2 e- iwt dw . (33) 
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B. Asymptotic behaviour at timelike infinity 



First, we consider the asymptotic behaviour of the fields at timelike infinity i + . As was 
explained, the late-time behaviour of the field should follow from the /o?/>frequency contri- 
bution to the Green's function. Actually, it is easy to verify that the effective contribution 
to the integral in Eq. ( |3"3"D should come from \w\=0(^). Thus, in order to obtain the asymp- 
totic behaviour of the field at timelike infinity (where y, y' <C t), we may use the |iu|r <C 1 
limit of $?i(r,w). Using Eq. 13.5.5 from \TM one finds 



\^i(r, w) ~ Ar 



1+1 



(34) 



Substituting this in Eq. ([33|) we obtain 
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Performing the integration in Eq. (|35|) , one finds 
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C. Asymptotic behaviour at future null infinity 

Next, we go on to consider the behaviour of the fields at future null infinity scri + . It is 
easy to verify that for this case the effective frequencies contributing to the integral in Eq. 
(|33|) are of order O(-). Thus, for y — y' <C t <C 2y — y' one may use the \vo\y' <C 1 asymptotic 
limit for ^/i(y',w) and the \w\y ^> 1 (Imw < 0) asymptotic limit of ^i(y,w). Thus, 

Ay n+ \ (37) 

and 
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where we have used Eqs. 13.5.5 and 13.5.1 of fUl , respectively. Integrating Eq. (|33"D with 

the aid of Eqs. ([37]) and (|38|) one finds (for t> 3> -u) 
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D. Asymptotic behaviour along the black- hole outer horizon 

Finally, we consider the behaviour of the fields at the black- hole outer- horizon r + . While 
Eqs. (^) and (E3j) are (approximated) solutions to the wave-equation ( |14D in the r 3> M, \a\ 
case, they do not represent the solution near the horizon. As y — > — oo the wave-equation 
(|T4|) can be approximated by the equation || 
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Thus, we take 



® 1 (y,w) = C(w)A- s/2 e 



s /I -iky 



(41) 



and we use Eq. (|37|) for ^i(y', w). In order to match the y <C — M solution with the 
y ^> M solution we assume that the two solutions have the same temporal dependence 
(this assumption has been proven to be very successful for neutral |19| and charged || 
perturbations on a spherically symmetric backgrounds). In other words we take C(w) to be 
w-independent. In this case one should replace the roles of y' and y in Eq. (|33|). Using Eq. 
(|33|), we obtain 
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where Tn is a constant. 
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E. Coupling of different multipoles 



The time-evolution of a wave-field described by Eq. (ffl) is given by Eq. (^). The 
coefficients Bi(r,9) and B 2 (r,9) appearing in Eq. (^) depend explicitly on the angular 
variable 9 through the rotation of the black-hole (no such dependence exist in the a = 
case). This angular dependence of the coefficients leads to an active interaction between 
different multipoles (characterized by different values of /). This coupling between different 
multipoles is physically caused by the rotational dragging of reference frames (due to the 
rotation of the black- hole). Thus, even if the initial data is characterized by a pure multipole I 
(i.e., by a certain spin-weighted spherical harmonic s Yj) other multipoles would be generated 
dynamically during the evolution. In other words, the late-time evolution of a spin-s field is 
dominated by the lowest allowed multipole, i.e., by the I — \s\ multipole, regardless of the 
angular dependence of the initial-data (provided that the initial-data contains modes with 
\ m \ < This phenomena of coupling between different multipoles has been observed in 
numerical solutions of the Teukolsky equation |5]|| (a similar phenomena is known in the 
case of rotating stars f20fl). 



V. SUMMARY 

We have studied the asymptotic late-time evolution of realistic rotating collapse. This 
was done by considering the asymptotic late-time solutions of Teukolsky's master equation 
which governs the evolution of gravitational, electromagnetic, neutrino and scalar (massless) 
perturbations fields on Kerr spacetimes. 

Following the no-hair conjecture for rotating black-holes we have shown that the asymp- 
totic solutions develop inverse power-law tails at timelike infinity, at null infinity and along 
the black-hole outer horizon (where the power-law behaviour is multiplied by an oscillatory 
term, caused by the dragging of reference frames at the event horizon). 

The damping exponents, describing the fall-off of the asymptotic solutions of Teukolsky's 
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equation at timelike infinity and along the black-hole outer horizon are independent of the 
spin parameter of the field. However, we have shown that the damping exponents at future 
null infinity are spin- dependent. Moreover, the asymptotic late-time solutions of Teukolsky's 
equation at all the three asymptotic regions have a spatial dependence on the spin parameter 
of the field. 

The damping exponents are independent of the rotation parameter a of the black-hole 
(or star). However, the rotational dragging of reference frames, caused by the rotation of the 
black-hole has two important effects (which are not found in the non rotating case) on the 
asymptotic late-time evolution of rotating collapse: The power-law tail along the black-hole 
outer horizon is multiplied by an oscillatory term. In addition, the rotation leads to an 
active coupling of different multipoles. Hence, the late-time evolution of a spin-s field is 
dominated by the lowest allowed multipole, i.e., by the I — \s\ multipole. 
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